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Abstract 

A heuristic model is given for anisotropic magnetohydrodynamics (MHD) turbulence in the 
presence of a uniform external magnetic field Bo en. The model is valid for both moderate and 
strong Bq and is able to describe both the strong and weak wave turbulence regimes as well as the 
transition between them. The main ingredient of the model is the assumption of constant ratio at all 
scales between the linear wave period and the nonlinear turnover timescale. Contrary to the model 
of critical balance introduced by Goldreich and Sridhar [P. Goldreich and S. Sridhar, ApJ 438, 763 
(1995)], it is not assumed in addition that this ratio be equal to unity at all scales which allows us 
to use the Iroshnikov-Kraichnan phenomenology. It is then possible to recover the widely observed 
anisotropic scaling law k\\ oc between parallel and perpendicular wavenumbers (with reference 
to Boe\\) and to obtain the universal prediction, 3a + 2(3 = 7, for the total energy spectrum 
E(k±,k\\) ~ In particular, with such a prediction the weak Alfven wave turbulence 

constant-flux solution is recovered and, for the first time, a possible explanation to its precursor 
found numerically by Galtier et al. [S. Galtier et al., J. Plasma Phys. 63, 447 (2000)] is given. 

PACS numbers: 47.27.Jv, 47.65. +a, 52.30.cv, 95.30.Qd 
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I. INTRODUCTION 



Turbulent flows are often studied under the assumptions of homogeneity and isotropy (see, 
for a review, Ref. 

ELQ) 

. Such assumptions are convenient for theoretical studies but are 
not always justified physically. For example, it is well known that stratification or rotation 
applied to neutral flows lead to anisotropic turbulence (see e.g. Ref. HQ]). Isotropy is 
even more difficult to justify in astrophysics where a magnetic field is almost always present 
at the largest scale of the system. The magnetohydrodynamics (MHD) approximation has 
proved to be quite successful in the study of a variety of space plasmas. During the last 
quarter of century many studies have been devoted to incompressible MHD turbulence in 



the presence o 
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a uniform external magnetic field B en (see Ref. 
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One of the most clearly established results is that the presence 
of B leads to a bi-dimensionalization of an initial isotropic energy spectrum: the turbulent 
cascade transfers energy preferentially to perpendicular wavenumbers, i.e. in the direction 
transverse to Bo ey . 

Constant-flux spectra are known to occur in many instances in turbulent flows, the best 
example of which being the Kolmogorov energy spectrum following a E(k) oc k~ 5 ^ 3 law 
for three-dimensional Navier-Stokes turbulence (see e.g. Ref. 2^). Power law spectra are 
also measured in turbulent MHD flows but the value of the scaling index is still hardly 
discussed in the community. The first prediction in MHD was given independently by 
Iroshnikov and Kraichnan in Ref. 21, 3] (hereafter IK). They argued that the destruction 
of phase coherence by Alfven waves traveling in opposite directions along the local large 
scale magnetic field introduces a new timescale and a slowing down of energy transfer to 
small scales. Assuming isotropy, the dimensional analysis for three wave interactions leads 
to a E(k) oc A; -3 / 2 spectrum for the total energy (see also Ref. Many direct numerical 

simulations of strong turbulence in isotropic (B = 0) MHD have been made during the 
last years (see e.g. Ref. 

but a definitive conclusion about the spectrum index 
is still not achieved mainly because the Kolmogorov and IK predictions are very close; 
furthermore, such scaling laws may be slightly altered by intermittency effects, and the 
numerical resolution is barely sufficient to determine such spectral indices. Goldreich and 
Shridar in Ref. proposed in 1995 a heuristic model of strong turbulence for anisotropic 
(moderate B ) MHD where the distinction between the perpendicular (k±) and parallel (k\\) 
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wavenumbers is made. This model is based on a critical balance between linear wave periods 
ta an d nonlinear turnover timescales tnl for which the equality ta = tnl is assumed to 
hold at all scales in the inertial range. (Only the symmetric case, for which the Alfven 
waves traveling in opposite directions carry equal energy fluxes, is considered here and in 
the remainder of the paper.) Then the one dimensional perpendicular power spectrum for 

r in 

the total energy scales as E(k±) oc k ± whereas the parallel and perpendicular spatial 
sizes of eddies are correlated according to the scaling law fcii oc k 2 / 3 . The latter prediction 

) and seems well observed in recent direct numerical 

f. 



simulations for moderate Bo (see Ref. 



121). 



As is well known, in the limit of large B MHD turbulence becomes strongly anisotropic 
and mainly weak in the sense that, in Fourier space, the domain of applicability of strong 
turbulence is confined to a region localized around the k\\ = plane. The formalism of weak 
Alfven wave turbulence developed by Galtier et al. in Ref. [lfl Q is well adapted to this 
situation. It leads to the so-called wave kinetic equations for which the exact power law 
solution for the total energy is E(k±,k\\) oc k'] 2 f{k\\). The function f{k\\) is undetermined 
because of the dynamical decoupling of parallel planes in Fourier space (this is the signature 
of the absence of energy transfer along the B direction; it thus represents a shadow of the 
initial conditions). Numerical simulations of the wave kinetic equations show clearly such a 
constant-flux spectrum but it also reveals the existence of a transient spectrum during the 
front propagation towards small scales with a steeper power-law in /c7 7 ^ 3 , the dynamics of 



which is not yet clarified (see e.g. Ref. m, 

l27j|). The discovery of such transient spectra 



in wave turbulence and the possible existence of a family of solutions that are not caught 
by the usual technique of conformal transform (given e.g. in Ref. j^sj]) constitute a new 
exciting 



Re 
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;opic of research where some progress is currently being made (see, for example, 
3Cj). When using a shell model of strong turbulence, it has also been found in Ref. 
that, when considering the decay of energy in time as a power law, initial transients occur 
that also follow power-laws to leading order and that precede the final power-law decay of 
the energy; the origin of such transients in time may well be a transient in the Fourier energy 
spectrum preceding the establishment of a Kolmogorov-like spectrum, although this point 
has not been documented yet. 

In this paper, we propose a heuristic model that describes anisotropic MHD flows for 
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the regimes of strong (moderate Bq) and weak wave (strong Bq) turbulence as well as the 
transition between them. As a result of our analysis, a family of solutions is found for 
the anisotropic total energy spectrum from which the transient spectrum described above 
is a particular solution. We show that the model supports the same anisotropic scaling 
law between the parallel and perpendicular wavenumbers as the one found in the context 
of critical balance but it is more general in the sense that here we do not impose equality 
between linear wave periods and nonlinear turnover timescales which allows us to use the 
IK phenomenology. We finally propose to extend the anisotropic model to two other types 
of fluids. 

II. ANISOTROPIC MHD MODEL 

The presence of the mean magnetic field Bq leads to anisotropy and thus to different 
variations with wavenumber directions, leading us to distinguish between kj_ and kn. More 
precisely we will assume that k± ^> k», i.e. that under the external agent (here, Bo), the 
turbulent flow develops principally in the direction perpendicular to this agent. We will 
consider the symmetric case for which, in particular, the r.m.s. value at any scale in the 
inertial range of the fluctuating velocity field v and magnetic field b have the same order 
of magnitude; note that b is taken in velocity units since we are restricting the analysis to 
the incompressible MHD case (the density is constant and can be assumed to be equal to 
unity) . 

In the classical Kolmogorov phenomenology (hereafter, K41), the fluctuations are dis- 
tributed isotropically and there is only one timescale, the nonlinear time or eddy turnover 
time tnl, which is also the transfer time of the energy to small scales within the system, 
i.e. r tr = tnl- The rate of energy transfer per unit mass writes Skai ~ E/r tr , where E 
is the total (kinetic plus magnetic) energy at a given scale. In the very same spirit, we 
can develop the IK phenomenology for anisotropic MHD turbulence using explicitly the fact 
that the eddy turnover time is t^l ~ (vk)~ l ~ (f&±) _1 and the Alfven wave period is 
ta ~ l/(k\\Bo). We assume that these two timescales are not equal which allows us to use 
the IK phenomenology. The rate of energy transfer per unit mass now writes 

f E v2 m 

tlK a ~ ~ , (J- J 

T tr T tr 
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where the transfer time will be assumed to scale as 



Ttr = TNL , {*) 



as is known to be the case for three- wave interaction processes (see e.g. Ref. [la, |2l|, |22|). 
The subscript "a" in IK a stands for the anisotropic version of IK. Including relation (J2J) into 
(JTJ, one obtains 

We now define the anisotropic energy spectrum, assuming self-similarity in both (_L, ||) 
directions: 

E(k±, k\\) ~ kj 01 , (4) 

where a and (3 are unknown. It is a 2D energy spectrum. In other words, the total en- 
ergy of the system is recovered by directly integrating the spectrum along kj_ and k\\, i.e. 
E sys — J J E(kj_, fc||) dk±dk\\ (see e.g. Ref. [r|, for a rigorous definition). It is important to 
note at this stage a difference with the classical phenomenology where the energy spectrum is 
not introduced with unknown indices but rather deduced by the analysis. The phenomenol- 



ogy proposed to describe M 
approaches in e.g. Ref. 0, 



ID turbulence is, of course, not unique; one can find recent 



. As is well known, the presence of an external magnetic 
field B leads to a reduction of nonlinear transfers along its direction which, in turns, leads 
to a reduction in size in Fourier space of the inertial range. Numerical simulations, for the 
most par t at moderate resolutions, and/or using the Reduced MHD approximation (see e.g. 
Ref. pj]) show indeed that a strong Bo leads to an absence of inertial range in k\\ with 

a spectrum mainly exponential (see e.g. Ref. o,y,y,Q) 

as it is usually observed in 

turbulence at large wavenumbers. Since heuristic models deal with power laws, they are 
mainly able to describe what happens inside the inertial range, without saying much about 
the strength of the nonlinear transfers and thus about the size of the inertial range. We 
are assuming, as noted before (see equation (4)) that power-law develops in k\\ as well, i.e. 
that the Reynolds number is extremely large, as found in astrophysical flows; in particular, 
it may need to be substantially higher than in the fully isotropic case since transfer in the 
fcii direction is strongly impeded in the presence of a uniform magnetic field. Therefore, 
relation (J1J) introduced above gives a description of the possible inertial ranges in both k± 
and h\ directions; additionally, it offers an opportunity to describe continuously the phase 
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transition between the strong and weak turbulence regimes. We introduce relation (j3J) into 
(jUJ), and by noticing that v 2 ~ E ~ E(k±, fcn) k±ku, we obtain: 

B E IKa ~k*- 2a k 1 - 2 ^ (5) 

which can also be rewritten in terms of a (k±, kn) relationship as: 



k^iBoSjK^k 2 /- 1 . (6) 

This is our first anisotropic relation. In order to proceed further, we are now seeking a 
second relationship between the two scaling indices a and (3. One way to obtain a unique 
relation between them is to use the assumption of constant ratio between t^l and ta at all 
scales in the inertial range. In other words, we will assume that 

X = ^ (?) 

TNL 

is a constant at all scales but not necessarily equal to one as it is in the critical balance model 
in Ref. [ill ]. Condition ((7|) can be understood as a formal balance between the linear and 
nonlinear terms in the ideal MHD equations. It is important to note that the precise value of 
X is not important in itself since it is just a numerical factor that does not change the power 
law scaling. However if x is smaller than one we are allowed to use the IK phenomenology 
instead of the K41 one for x = 1. A ratio of one seems to be very restrictive and does not 
correspond to some of the results stemming from direct numerical simulations where x can 
be smaller than unity, as observed for example in Ref. a two-dimensional geometrical 
case where local anisotropy is possible, or in solar wind in situ measurements (see Ref. j^) 
where x seems to be smaller than unity. It is also in apparent contradiction with wave 
turbulence theory for which we have \ <C 1. We will see that the assumption of constant 
ratio at all scales in the inertial range is sufficient to achieve a unified model that describes 
both strong and weak MHD turbulence. With the definition of the timescales given above, 
we obtain 

X~—. (8) 

X B k\\ 1 ' 



Including expression (JHJ) into @ gives: 
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If we use explicitly the fact that both y and the rates of energy transfer per unit mass do 



not depend on the scale (see e.g. Ref. [20 

fci 



21 



we finally obtain 



k 



2/3 



Bq 



(10) 



This leads us to the first main conclusion of this paper: the same scaling in wavenumbers 
as the critical balance model is obtained when \ is not equal to unity, i.e. when the IK 
phenomenology is used. It thus seems to be a general law for incompressible MHD turbu- 
lence, when either the Kolmogorov phenomenology for energy transfer or the IK a three-wave 
formulation of energy transfer is utilized, as long as a critical balance (generalized to any 
ratio differing from unity) is assumed between stretching by velocity gradients and wave 
motions in the presence of a uniform field Bq. Condition ()10|) may be seen as a path in 
k±-k\\ Fourier space followed naturally by the dynamics for different values of the (imposed) 
uniform magnetic field. An illustration is given in Figure for two different values of 
Bq. It is interesting to note that condition (jlOj) is not incompatible with the weak wave 
turbulence prediction since in the infinite limit of Bq the curve identifies with the k±-a.xis 
and therefore no transfers are possible along the parallel direction. 



k//i 



2/3 




FIG. 1: Path, in Fourier space, followed by the MHD turbulence dynamics for two different values 
of the external magnetic field Bq. A higher value of Bq leads to weakening of parallel transfers 
with a curve closer to the k±-axis. When Bq is infinite the path identifies with the fcj_-axis and no 
transfers are possible along the parallel direction. Note that the curves are not extended down to 
the origin since the anisotropic law is only valid in the inertial range. 



Another important and somewhat unexpected consequence of this analysis has to deal 
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with the scaling of the energy spectrum. In the context of critical balance, it is claimed 
that the energy spectrum derived with an anisotropic Kolmogorov phenomenology scales 
like E(k±) ~ k^ 5 ^ 3 . Here we show that it is in fact possible to find a multitude of spectral 
indices a and j3 that satisfy the assumption of constant ratio x at all scales in the inertial 
range. Such indices follow in fact a general linear relationship. To find this relation we need 
to equalize the power-law behaviors found in relations © and (|10|) between the parallel and 
perpendicular distribution of modal energy in Fourier space. If we use explicitly the fact 

hat the rates of energy transfer per unit mass do not depend on the scale (see e.g. Ref. 

2Q, 21, we finally obtain 

3a + 2/5 = 7. (11) 

Relationship (|TT|) is general and can be used for strong turbulence as well as for weak wave 
turbulence. Note that the relationship ()11|) is compatible with the anisotropic Kolmogorov 
spectrum (with a = 5/3 and /3 = 1) as advocated by Goldreich and Sridhar in Ref. [llj], and 
with the anisotropic IK-like spectrum corresponding to three-wave interactions for weak 
MHD turbulence {a = 2 and (3 = 1/2) (see Ref. Q, 40, 41|). We conjecture that it 



is compatible with the physics of the transitional regime between weak and strong MHD 
turbulence as well. In other words the law (fTTj) we just derived shows that the 2/3 scaling 
between a dependence of k±_ and fcii is not a unique signature of the Kolmogorov spectrum 
but rather a signature of the rate at which energy is transfered to small scales where it is 
dissipated, i.e. a trace of the decay of energy, and as such an unavoidable dimensional law 
assuming equation (@J) holds and that \ is independent of wavenumbers in the inertial range. 

III. DISCUSSION 

The preceding remark may be linked to the heretofore unexplained following finding: 



using direct numerical simulations in two space dimensions, it was shown in Ref. 



421 that 



scales (as expressed in 



44|) have a self-similar 



the structure functions of order p based on the energy flux to smah 
terms of the exact laws for MHD turbulence derived in Ref. [ 
scaling in the inertial range which is compatible with the scaling for the velocity structure 
functions in fluid turbulence, whereas the structure functions of the basic fields (velocity 
and magnetic fields) are more intermittent insofar as they depart more significantly from a 
linear scaling with p. What this paper shows is that both models (K41 and IK) can be seen 
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in a unifying way, illustrated by the law (fTTj) . Note that in the case of the advection of a 
passive tracer such as temperature, it also seems that the scaling of the structure functions 
based on the flux of the tracer is close to the fluid scaling laws (see Ref. whereas 
the tracer itself is well-known to be strongly intermittent. 

More surprisingly perhaps, the choice of a = 7/3 and j3 = also satisfies the law 
derived in this paper. This k ± 7 ^ 3 spectrum was found in Ref. Q| as a precursor to the 
constant flux solution of the Alfven wave kinetic equations which establishes itself later in 
time; indeed, in that paper, the 7/3 spectrum is found numerically for the case kn = 0, 
i.e. without any dependence in k\\ (corresponding to /3 — 0). The simple model proposed 
here thus sheds some light, albeit heuristically, on two intriguing facts that have emerged 
recently concerning weak wave turbulence: (i) the fact that, preceding in time the constant 
flux solution, a power-law spectrum (called the precursor) establishes itself the origin of 



which was unrelated to anything known about turbulence spectra (see Ref. lg, |29j) until 
now, and (ii) the fact that in some cases (see e.g. Ref. 3Cj) in the weak wave turbulence 
regime, a wealth of power-law solutions in the (oc,/3) plane can be found numerically as 
stationary solutions to the wave kinetic equations, although it is not clear whether such 
solutions are attractive, nor whether they are stable. The link between the dimensional 
argument given here which allows to recover the a = 7/3 precursor spectrum, and the 
self-similar argument given in Ref. Q|, which is compatible with the a = 7/3 spectrum, 
remains to be clarified. Indeed, in Ref. one recovers the numerically observed law of 
self-similar decay of the energy spectrum, assuming that the energy spectrum scales in the 

— 7/3 

precursor phase as k ± , whereas in this paper we give a heuristic justification to that same 
—7/3 law in terms of dimensional arguments compatible with three- wave interactions and 
with the assumption that there is no k\\ dependency in the precursor solution. We show 
furthermore that this law stems from the same type of unified approach that also gives the 
other known spectra for MHD turbulence. In order to ascertain the validity of the model 
derived here, numerical computations with the full three-dimensional MHD equations can 
be envisaged but being able to distinguish between different power laws which are close in 
their respective spectral indices may prove difficult. This point is currently investigated and 
will be presented elsewhere. 

The argumentation delineated here could be used for other types of wave turbulence, 
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like in the case of whistler waves (see Ref. Q), inertial waves (see Ref. Qj) or gravity 
waves (see Ref. F° r whistler waves that can be encountered in Hall MHD or electron 

MHD, the characteristic wave period scales as Tw ~ (see e.g. Ref. One finds, 

bliowing the same analysis as presented before in this paper, that ku ~ k 1 / 3 (see also Ref. 

5o| ) : hence, the prediction for scaling exponents as defined in equation (@J) for the whistler 
turbulence energy spectrum becomes: 3a + (3 = 8. (Note that here, the eddy turnover time 
is based on the magnetic field, i.e. t^l — [v = since v oc V x b.) The known 

constant-flux solution, a = 5/2 and (3 = 1/2, to the wave kinetic equations derived in Ref. 

is recovered again as well as the strong turbulence prediction, a = 7/3 and (3 — 1 (see 
e.g. Ref. Q)- 

When considering inertial waves in rotating turbulence for a Navier-Stokes fluid, the wave 
period scales as 77 ~ k±/ku. The same analysis as before leads to the following relationship 
between scaling exponents: 3a + 5(3 = 10. The known constant-flux solution to the wave 
kinetic equations written in Ref. 4^ for wave turbulence corresponds to a = 5/2 and (3=1/2, 
which does fulfill the above relationship and we recover the Kolmogorov prediction for strong 
turbulence as well, viz. a = 5/3 and (3 = 1 for a weak rotation rate (see the comment in 
Ref. 0)- 

This prediction again could be tested using direct numerical simulations, and the 
possible existence of precursors could also be studied and checked against the compatibility 
with the above relationship between the scaling exponents a and (3 in the rotating case. 
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